Observation of a warped helical spin-texture in Bi 2 Se3 from circular dichroism 
angle-resolved photoemission spectroscopy 
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A differential coupling of topological surface states to left- versus right- circularly polarized light is 
the basis of many opto-spintronics applications of topological insulators. Here we report direct evi- 
dence of circular dichroism from the surface states of Bi2Se3 using a laser-based time-of- flight angle- 
resolved photoemission spectroscopy. By employing a novel sample rotational analysis, we resolve 
unusual modulations in the circular dichroism photoemission pattern as a function of both energy 
and momentum, which perfectly mimic the predicted but hitherto un-observed three-dimensional 
warped spin-texture of the surface states. By developing a microscopic theory of photoemission 
from topological surface states, we show that this correlation is a natural consequence of spin-orbit 
coupling. These results suggest that our technique may be a powerful probe of the spin-texture of 
spin-orbit coupled materials in general. 



Three-dimensional topological insulators [lj— are an 
intensely researched phase of matter owing to their 
unique spin- helical metallic surfaces, where the spin di- 
rection is locked to the wavevector and winds by 2tt 
around the Fermi surface (j-[6|. Spin-dependent absorp- 
tion of circularly polarized light has recently been pre- 
dicted to generate highly spin-polarized surface electrical 
currents whose direction can be switched by the light he- 
licity However, a demonstration of spin-dependent 

differential absorption of left- versus right-circularly po- 
larized light [circular dichroism (CD)] from the surface 
states (SS) has evaded conventional probes such as trans- 
port and optics because of the combined need for energy- 
momentum resolution and surface sensitivity. 

In this Letter, we map the variation of the CD in the 
photoemission intensity jlol-[l5| over the full surface band 
structure of the topological insulator Bi2Se3 [HI, UJj us- 
ing a laser-based time-of- flight angle-resolved photoemis- 
sion spectroscopy (TOF-ARPES) technique. Our results 
show a strong CD from the SS in agreement with a recent 
work on Cua3i2Se3 [18j. However, owing to our ultra- 
high polarization purity and our novel sample rotational 
analysis, we uniquely resolve fine modulations in the CD 
pattern as a function of energy and momentum that per- 
fectly follow the predicted 19, 20|, but so far un-detected 
three-dimensional spin-texture of Bi2Se3. By devel- 
oping a microscopic theory of photoemission from topo- 
logical SS, we show that such modulations naturally arise 
from the coupling of the circularly polarized light to the 
surface spin- texture via spin-orbit interactions. Excellent 
agreement between our results and theory suggests that 
our technique is an efficient and sensitive measure of all 
three components of the surface spin-texture. 

The experimental setup is shown in Fig. 1(a). Single 
crystals of Bi2Se3 are cleaved along their (111) surface at 
room temperature under ultra-high- vacuum (< 6 x 10 -11 
Torr). Femtosecond laser pulses from a TiiSapph ampli- 
fier are frequency quadrupled to 6.2eV and then circu- 
larly polarized to > 99% purity with a quarter wave- 



plate. Strain induced birefringence from the vacuum 
windows and dichroism from all optical elements follow- 
ing the quarter wave-plate, which can ruin the purity of 
circular polarization, were carefully compensated in our 
experiment [21[ . Conventional hemispherical ARPES an- 
alyzers rely on the deflection of electron paths along one 
spatial dimension to measure electronic dispersion along 
one direction in momentum space. In contrast, our TOF- 
ARPES analyzer is developed for pulsed photon sources 
and resolves electron energies through their flight time 
from sample to detector, which enables measuring dis- 
persion spectra I(E,k x ,k y ) simultaneously over two di- 
mensional momentum space. Complete spectra of Bi2Se3 
(surface Dirac cone plus bulk conduction band) obtained 
with TOF-ARPES is shown in Fig. 1(b). This ability to 
measure the entire bandstructure simultaneously elimi- 
nates experimental geometry induced matrix element ef- 
fects [22] that may obscure the intrinsic CD pattern and 
allows for direct comparison of CD-ARPES intensities 
across different phase space points. 

To investigate whether CD is exhibited by Bi2Se3 and 
if it has bulk and/or surface state origin, we measure 
spectra like in Fig. 1(b) using left- and right-circularly 
polarized light and then take their difference to obtain 
the CD spectrum AI(E, k x ,k y ). As shown in Figs. 1(c)- 
(d), there is clear CD from the surface Dirac cone of 
Bi2Se3 and not from the bulk conduction band. The 
CD pattern exhibits perfectly odd reflection symmetry 
about the photon scattering plane when it is parallel to 
the YM direction [Fig. 1(e)], which is expected from 
crystal mirror symmetry. However this odd reflection 
symmetry is lifted when the scattering plane is along a 
non mirror-symmetric direction such as YK [Fig. 1(f)], 
which is evidence that the CD pattern is intrinsic to the 
sample. 

Having shown that CD is intrinsic to the sample and 
present for the SS, we first focus on its fine features in 
the low energy region [±0.1eV relative to the Dirac point 
(DP)] where the Fermi surface is circular. Fig. 1(d) 
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FIG. 1. (a) Schematic of a time-of- flight based angle- resolved 
photoemission spectrometer. PSD: position-sensitive detec- 
tor, (b) A typical iso-intensity surface in (E,k x ,k y ) space 
from Bi2Se3 collected simultaneously using linearly polarized 
photons, (c) Difference of TOF-ARPES data measured using 
right- and left- circularly polarized light, hv denotes the inci- 
dent photon direction, (d) E — k cut through the CD data 
volume in (c) along the angle denoted by the green dash in (e). 
(e-h) Constant energy slices through (c) at the Fermi level and 
at O.leV above the DP (insets) for sample rotation angles of 
= 0°, 30° , 60° and 90° respectively. Green hexagons repre- 
sent the Brillouin zone of Bi2Se3(lll), red lines are the mirror 
planes and the arrow denotes the photon incident direction. 



shows that at a constant energy the sign of CD is re- 
versed for states of opposite momentum and is reversed 
again across the DP. Furthermore, the CD pattern is in- 
variant under rotation ((/>) of the crystal mirror plane 
relative to the scattering plane [Figs. 1(e)- (h) insets]. 
These behaviors match the known in-plane spin texture 
in the low energy region, which are tangential to the ro- 
tationally isotropic constant energy contours 0, 0] and 
reverse chirality across the DP. 

The CD patterns change drastically away from the 
low energy region (\E\ > O.leV), where the Fermi sur- 
face evolves from being circular to hexagonal in shape 
Til . 23[. Here we observe new modulations in the magni- 
tude and sign of CD as a function of the angle 6 around 
the constant energy contours [Figs.l(e)-(h)]. In partic- 
ular, extrema along the k y axis (at = 90° and 270°) 
at low energy [Figs. 1(e)- (h) insets] become nodes in the 



high energy <ft = 0° spectrum [Fig. 1(e)]. Furthermore, 
in stark contrast to the low energy region, CD patterns 
at high energy undergo dramatic changes as <p is varied. 
For instance, the nodes along k y in the Fig. 1(e) spec- 
trum become extrema again in the <j) = 60° spectrum 
[Fig. 1(g)] and only repeats itself under <p — 120° rota- 
tions, consistent with the symmetry of the underlying lat- 
tice. This is reminiscent of the predicted dependence of 
the out-of-plane component of spin (S z ) at high energies 



19, 



201 ]. which goes from minimum to maximum under 
<p = 60° rotation. Such behavior has been observed us- 
ing Mott polarimetry only in Bi2Te3 [5] but not in Bi2Se3 
[6[ owing to the much larger value of (S z ) in the former. 
Given that the overall CD patterns observed closely fol- 
low the predicted spin-texture of the SS, we investigate 
their possible connection by developing a standard model 
of photoemission from topological SS. 

The microscopic Hamiltonian for a system with spin- 
orbit coupling is given by: 



H 



p2 



r(P x W) 



(1) 



where P is momentum operator, V(f) is the crystal po- 
tential, and s is the electron spin. Coupling to an elec- 
tromagnetic field is obtained via P — » P — eA, where A 
is the photon vector potential, such that to first order in 
A: 



H(A) = H -V A 



(2) 



where V = 4r ^2 c 2 (VV xs). The photoemission ma- 

trix element between the initial and final states is given 

by 

M(k,f) = (f j:\f-A\k) (3) 

where A = J dtA(t)e lUJt is the Fourier transform of 
A and \f^) is the bulk final state. The initial state 

|fc) = u^\(j)\) + v^\(j) % _) is a linear combination of two- 
fold degenerate pseudospin states \(j) l ±) at the DP that 
are eigenstates of total angular momentum (orbital plus 
spin), which is widely used in standard k • p descrip- 
tions of topological SS @, E3, 0, E3- Because of strong 
spin-orbit coupling, only pseudospin is a good quan- 
tum number [19[. The coefficients and deter- 
mine the expectation value of three pseudospin compo- 
nents: (S x ) % = h(ulv^ + vtu^, (S y ) n = h(vtu^-iulv^, 
{S z }^ = h(\u^\ 2 — \v^\ 2 ). Importantly, because spin 
is directly proportional to pseudospin [19], we refer to 
the two interchangeably. For circularly polarized light 
incident onto the surface with wavevector in the xz 
plane, A(t) = (A x sinc^t, A y cosc^t, A z sin uot) and A = 
(— \A X , A yi — \A Z ). Straight-forward application of time- 
reversal and crystal symmetries [2l[ yield the following 
expression for the photoemission transition rate: 

I(k) = a 2 (\A X \ 2 + \A y \ 2 ) + b 2 \A z \ 2 + a 2 lm (A x A* y ) (S Z ) R 
+ 2ablm[A* x A z (S y )f: - A* y A z (S x ) R } (4) 
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FIG. 2. (a) x, (b) y and (c) z components of the spin-texture 
over the complete surface states obtained by summing or sub- 
tracting AI data volumes at different 4> (see text). Color maps 
were created by integrating the data radially in k-space over a 
±0.015A _1 window about the surface state contours at each 
energy. Constant energy cuts at the Fermi level for each spin 
component are shown directly below in (d) through (f). 



where a and b are bandstructure dependent complex con- 
stants and Im refers to the imaginary part. Circular 
dichroism is obtained by taking the difference of Eq.(j4]) 
with opposite photon helicity (A y — )► — A y ): 

AI = a 2 Im(A x A;)(S z )t - Aablm(A z A* y )(S x ) % (5) 

Although the CD is a linear combination of (S x ) and 
(S z ), these can be disentangled by applying symmetry 
properties of Bi2Se3(lll) as follows. Time-reversal sym- 
metry and three-fold rotational symmetry together dic- 
tate that (S z ) flips sign upon a A(j) = 60° rotation while 
(S x ) stays unchanged. Taking the difference (sum) of 
CD patterns Acj) = 60° apart isolates (S z ) ((S x )). The 
(S y ) component is trivially obtained by performing the 
procedure for (S x ) under a 90° sample rotation. Fig. 
2 shows all three spin components measured using this 
method. Because calculating a and b is beyond the scope 
of our work, we only plot the relative and not absolute 
magnitude of the spins. We immediately notice that all 
three components reverse sign across the DP as expected 
and that (S x ) and (S y ) are modulated with a similar pe- 
riodicity. On the other hand (S z ) exhibits a different 
periodicity at high energies, which vanishes in the low 
energy region. 

To quantitatively test the validity of our microscopic 
theory, we first compare our extracted spin components 
with the ideal spin texture in the low energy region that is 
experimentally well-known . Constant low energy slices 
through the spin maps [Figs. 2(a)-(c)] show that (S z ) ex- 
hibits negligible modulation whereas (S x ) and (S y ) fol- 
low a clear sin(0) dependence that reverses sign across 
the DP [Figs. 3(a)-(c)]. By projecting the planar spin 
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FIG. 3. The evolution of the (a) S x , (b) S y and (c) S z compo- 
nents around constant energy contours in the vicinity of the 
Dirac point, obtained by taking cuts through the data in Figs. 
2(b)-(d). Curves are fits to sin(0) and cos(#) varying functions 
for S x and S y respectively, (d) and (e) show the spin vectors 
(each normalized by its magnitude \S\) of all positive and neg- 
ative energy states within zbO.leV respectively projected onto 
the k x — k y plane. 



vectors of each state in the low energy region of the Dirac 
cone onto the x — y plane, we obtain a vector field that ex- 
actly matches an ideal helical spin-texture with opposite 
chirality for electrons and holes [Figs. 3(d) and (e)]. We 
note that while the relative orientation of spins is directly 
measured in our experiment, the absolute sense of chiral- 
ity was set by matching our data to Mott polarimetry 
results at a single energy- momentum point [6] , although 
it can also be independently obtained by calculating a 
and b. The excellent agreement between our measured 
low energy CD maps and theoretical [19| and Mott po- 
larimetry results [j-lSl strongly suggests that CD-ARPES 
is a sensitive measure of the topological SS spin-texture. 

Next we compare our CD spectra to our model in the 
high energy region where the spin-texture is predicted to 
depart from ideal planar helical form but is so far un- 
observed in Bi2Se3 using conventional Mott polarimetry 
based ARPES [6] . At high energies where the Fermi sur- 
face acquires a hexagonal shape, k ■ p theory predicts 
that the spin-texture should develop a finite out-of-plane 
component (S z ) that is modulated with a sin(3#) period- 
icity, with maximally upward and downward spin cant- 
ing along YK directions [19]. Fig. 4(a) shows a con- 
stant energy slice at 0.3eV through the spin map of Fig. 
2(c), where there is a clear sin(3#) dependence of (S z ) 
that is maximal along f K and zero along f M, in per- 
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Deformed Dirac spin-texture 

FIG. 4. Evolution of (a) S z and (b) S x around the Fermi 
surface. Red line in (a) is a fit to k • p theory. The blue and 
pink lines in (b) are fits to phenomenological models 21] that 
respectively do and do not account for an in-plane angular 
deviation (5) away from S JL k locking behavior (inset), (c) 
The energy dependence of the amplitude of S z modulation 
(green symbols), its fit (magenta line) and 5 (red symbols). 
Inset shows the spin- integrated Fermi surface together with a 
fit to k-p theory, (d) Spin-orientation around the hexagonally 
deformed Fermi surface constructed from CD-ARPES data. 
3D view shows spins plotted around Fermi contour and top- 
down view shows spins projected onto a segment of the unit 
circle in momentum space for ease of visualizing 5. 

feet accord with both theory [19] and experiment [H on 
Bi 2 Te 3 . To further test whether our CD-ARPES data is 
indeed sensitive to spin, we compare our measured en- 
ergy dependence of out-of-plane spin component along 
YK ((S Z )°(E)) with its unique functional form predicted 
by k • p theory [l9| : 



(S z )°(E) = l/^l + [k(E)(3}- 



(6) 



where j3 is the only free parameter. Eq.(|6]) provides an 
excellent fit for our data [Fig. 4(c)] and the fitted value 
of /3, which parameterizes the degree of hexagonal distor- 
tion of the Fermi surface, is quantitatively consistent with 
the parameters extracted from fitting the Fermi surface 
contour in Fig. 4(c) inset 21]. The accurate extraction 
of all three spin components over the entire Dirac cone 
from our CD-ARPES data [Fig. 4(d)] is a testament to 
the validity and sensitivity of our method. 

While our data produce all features predicted by k • p 
theory [l9|, the in-plane spin component at high en- 
ergies also exhibits unpredicted kinks in its 6 depen- 
dence [Fig. 4(b)]. One possible origin of these kinks 
is that the in-plane spin magnitude is modulated purely 
by the out-of-plane spin canting. However a fit to the 
data using this assumption [Fig. 4(b) pink line] yields 
kinks that are too weak and fails to reproduce the nodes 
at = 0° even when allowing for unrealistically large 
(~ 45°) canting angles [2l|. An alternative explanation 



is that the in-plane spin direction deviates from perpen- 
dicular spin-momentum locking [inset Fig. 4(b)] period- 
ically in 6. A fit assuming periodic angular deviations 
of amplitude S [Fig. 4(b) blue line] better reproduces 
the data. The in-plane component constructed from the 
data alone [Fig. 4(d)] is consistent with the latter in- 
terpretation and shows that the deviation is zero along 
YM and YK and is maximum in between the two. A 
plot of the amplitude of S versus energy [Fig. 4(c)] shows 
that it only develops beyond ~ 0.3eV from the DP. The 
spike exactly at the DP is consistent with spin degener- 
acy at the DP. This in-plane spin canting behavior has 
been shown by first-principles calculations as well as k ■ p 
expansions up to fifth order on Bi2Te3 [24[. However, 
although such effects have been observed in Rashba spin 
split metals [25j, they have never been observed in any 
topological insulators. 

In conclusion, we have directly observed circular 
dichroism from the surface states of Bi2Se3. By com- 
bining a TOF-ARPES technique with ultra-pure circular 
photon polarization and a novel sample rotation analysis, 
we resolve unusual modulations in the circular dichro- 
ism photoemission pattern as a function of both energy 
and momentum for the first time, which closely follow 
the predicted three-dimensional spin-texture. A direct 
connection between CD-ARPES and spin-orbit induced 
spin-textures is established through our microscopic the- 
ory of photoemission. Our results open the possibility 
to generate highly-polarized spin currents with circularly 
polarized light, which may be detected through transport 
or optical means [26]. The efficiency and high spin sen- 
sitivity of our technique suggest that CD-ARPES may 
be used as a vectorial spin mapping tool to detect small 
deviations from a tt Berry's phase in magnetically doped 
topological surface states, or to study spin-orbit coupled 
materials in general. 

Towards completion of this manuscript, we became 
aware of three related concurrent and independent works 
18 , 27( | . None, however, resolve the fine modulations in 
the CD-ARPES patterns that are integral to understand- 
ing the connection to spin, which is the main focus of this 
work. 
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